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O . Abstract 

cr 

5^ | Properties pertaining to thermodynamical local stability of Reissner-Nordstrom black holes sur- 

rounded by quintessence as well as adiabatic invariance, adiabatic charging and a generalized Smarr 
formula are discussed. Limits for the entropy, temperature and electric potential ensuring stability 
of canonical ensembles are determined by the classical thermodynamical and Poincare methods. By 
0^ . the latter approach we show that microcanonical ensembles (isolated black holes) are stable. Two 

geometrical approaches lead to determine the same states corresponding to second order phase tran- 



l/"") \ sitions. 
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1 Introduction 



X 

Quintessence has made the subject of many papers ranging from exact solutions ul to cosmological 
models []2] (and references therein) all fueled by observational data, which were achieved through various 
projects 0, and the discovery of the acceleration of the universe. These observation lead to believe that 
the accelerated expansion is due an exotic fluid with negative pressure making up the quintessence. 

The spherically symmetric static solutions derived in (TJ, with quintessence as source term, obey a 
special condition of additivity and linearity in the energy-momentum tensor. They do not exhaust the set 
of spherically symmetric static solutions and the quest for new solutions remains open. 

The solutions derived in [1] depend on four parameters: the mass and charge (M,q) of the Reissner- 
Nordstrom (RN) black hole and (c, co) where c > 0, which might be called the charge of quintessence, 
determines the energy density p q of quintessence and — 1 < CO < relates pressure to energy density 
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via the equation of state p q = (op q . A first classification of these solutions is based on their asymptotic 
behavor: 

— 1/3 < © < 1 : asymptotically flat solutions 

— I < (O < — 1/3: non-asymptotically flat solutions . 

Only the physical, geometrical and thermodynamical properties of the non-asymptotically flat solu- 
tion with CO = — 2/3 were investigated in [1J. One of the aims of this paper is to discuss those physical 
and thermodynamical properties of asymptotically flat solutions (— 1/3 < CO < 1) pertaining to thermo- 
dynamical stability (Sections |2] and [3]). We mainly discuss: horizon location, extremality conditions, 
relevant thermodynamic entities, generalized Smarr formula and first law of thermodynamics, adiabatic 
invariance and adiabatic charging of RN black holes with or without quintessence. Non-asymptotically 
flat solutions will make the subject of a subsequent work. 

Gravitating systems do not obey the linear rules for mass and entropy additions. Thermodynamical 
methods appealing to ensembles where mass or entropy is held constant (and used as a control parameter) 
do not apply to gravitating systems SH6]. An instance of that, the thermodynamic stability analysis or 
phase transitions of an isolated black hole (mass does not fluctuate) can't be carried out by either classical 
thermodynamical or geometrical methods which rely on the linearity hypothesis. 

Poincare Q developed a powerful method applicable to problems pertaining to equilibrium and con- 
ditions of stability. Originally the method, known as the turning point method (TPM), was applied to the 
uniform rotational motion of a homogeneous liquid to determine the cases of local equilibrium and the 
conditions of stability of such equilibrium. Then it was applied to different situations ||8l- ifTBl includ- 
ing problems non-tractable by classical thermodynamical or geometrical methods. It was generalized to 
many -parameter equilibrium families |[T4ll . 

It is always possible to subdivide the space of all equilibrium configurations into 1 -parameter sub- 
spaces. In each subspace, all points representing equilibrium states are related by varying one parameter 
called the control parameter; hence the name of linear series of equilibrium given to each subspace. The 
method of Poincare [7 ] consists in further subdividing each linear series of equilibrium into smaller sub- 
spaces labeled stable, less stable, unstable states of equilibria. The method employs the terminology 
of "degree of stability", increasing or decreasing depending on the number of negative modes of the 
associated Hessian matrix, and no notion of phase transition is employed. The choice of the Massieu 
function or the thermodynamic potential depends on the thermodynamic ensemble under consideration 
and in many cases there is no need to evaluate the eigenvalues of the Hessian matrix associated with the 
Massieu function or the thermodynamic potential (H- lfl2l . 

Section |4] is devoted to the determination of the conditions of stability of equilibrium configurations 
of RN black holes surrounded by quintessence. We apply both, but in reversed order, the TPM and, 
whenever applicable, the classical thermodynamic approach. We will be able to always reach the same 
conclusion regarding the stability conditions but using different terminologies. 

Geometrical methods lfT5l - |fl9l are geometric approaches which attach a measure of length to the 
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space of all equilibrium configurations, the metrics of which are built up from a Legendre-invariant 
thermodynamic potential and its first and second order partial derivatives with respect to a set of ex- 
tensive variables. It is well known that stability analysis results depend on the thermodynamic ensem- 
ble Bti6ll8l[T0l, despite this fact the results derived in |fT5l - lfT8 1 do not make reference to any ensemble. 
Thermodynamic ensembles were used for the first time in [19]. The geometric methods rely on the lin- 
earity hypothesis and so are not applicable to isolated self gravitating systems. In Section [5] we define, 
within the contexts of these geometric approaches, the canonical ensemble describing the thermodynamic 
of an RN black hole immersed in a heat bath then apply the methods to analyze its stability. We conclude 
in Section [6l 

2 Properties of RN black holes in presence of quintessence 

In 4-dimensional spacetime, a spherical symmetric RN black hole plunged into the field of a spherical 
symmetric quintessence has the metric (H 

ds 2 = f(r)At 2 - f-\r)Ar 2 -r 2 m 2 (2.1) 

witrQ 

2M a 2 2c 

/(r) = 1 __ + 72___, -KflXOandoO. (2.2) 
With this notation, the density of energy of quintessence is 

3(Oc 

Note that the integration of p q over a hyper-surface of constant t diverges as lim r ^oo r~ 3w . As we shall see 
later, this results in a similar divergence of the Komar mass of the black hole. This infinite background 
contribution due to the density of quintessence should be subtracted in an appropriate way so that the 
Komar mass of the hole remains finite. 

The physical properties of RN black holes described by (12.11 ) depend on the sign of 3(0 + 1. We shall 
consider the case where the asymptotic behavior of the hole is not altered by the presence of quintessence, 
that is the case where the black hole solution is still asymptotically flat. This correspondences to 3(0 + 1 > 
(— 1/3 < to < 0) with further constraints as shown below. 

In the case 3(0 + 1 = 0, the metric (12.11 ) is not always asymptotically flat. If we assume M > 
then 1) if 1 — 2c > (0 < c < 1 /2), the metric may be brought to the following form upon performing 
the coordinate and parameter transformations: t' = \/l — let, r' = rj \f \ — 2c, M' = M/(\ — 2c) 3 / 2 and 
q' = q/{\-2c): 

As 2 = f'dt' 2 - f'- l dr' 2 - (1 - 2c)r a &a 2 (2.4) 

where /' = 1 — 2M' jr 1 + q' 2 /r' 2 . This asymptotically flat metric (12.41 ) has a conical singularity in each 

'in the original derivation of d2.2t . c was taken negative flj. We have made the substitution — c — > 2c for simplicity and 
convenience. 
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plane 6 = do = constant. The deficit angle depends on do and is equal to 4nc in the plane 6 = n/2. 2) if 
1 — 2c < 0, the metric d2.il ) is no longer asymptotically flat. 

Let -l/3<ft)<0(l>3ft) + l>0). The horizons are defined by the condition f(r) = 0. Setting 
u=\/r, this implies 

\-2Mu + q 2 u 2 = 2cu 3w+1 . (2.5) 

The parabola y = 1 — 2Mu + q 2 u 2 has an absolute minimum value of (q 2 —M 2 )/q 2 at w m i n = M/q 2 . 
Thus, if q 2 < M 2 the graphs of y = 1 — 2Mu + q 2 u 2 and y = 2cu 3co+l always intersect, as shown in 
Figure Q] (a), at two points («/j,wi) such that ut, < u + and u\ > u_ with 1/m+ = r + = M + \JM 2 — q 2 
and 1/m_ = r_ = M — \JM 2 — q 2 . The event horizon is at r/, = l/uf, > r + . Here is another effect of 
quintessence: For fixed M and q, the radius of the event horizon increases leading to higher entropy with 
respect to that of an RN black hole. 



Y (a) 



Y (b) 



Y (c) 




Figure 1 : Plots of y = 1 - 2Mu + q 2 u 2 (dotted line) and y = 2c« 3<B+1 for - 1/3 < ££>< 0. (a) A black hole: q 2 <M 2 . The two 
curves intersect at two points labeled (b) A black hole: q 2 > M 2 and c > c c . The two curves intersect at two points 



labeled (uh,ui). (c) An extreme black hole: q > M and c = c c . The two curves intersect at the point uq given by 02. 9t . 

Even in the case q 2 > M 2 , the two curves y = 1 — 2Mm + g 2 w 2 and j = 2cu 3co+l may still meet at two 
points (uh,u\) as shown in Figure [T](b). Hence, the presence of quintessence makes it possible for the 
black hole to absorb more charge than the ordinary RN black hole before it becomes a naked singularity. 
For fixed ft), M and q 2 /M 2 > 1 there is a critical value c c of c below which the two curves do not intersect. 
c c is such that the two curves y = 1 — 2Mu + q 2 u 2 and y = 2cu 3co+l have a common tangent line at the 
unique point of intersection uo > 0: 



1 - 2Muo + q 2 uo 2 = 2c c uo 3(0+1 
-M + q 2 uo = c c (3ft) + l)wo 3f0 . 



(2.6) 
(2.7) 



Eliminating c c uo ia leads to 



(1 - 3co)q 2 uo 2 + 6Mft>M - (3ft) + 1) = 



(2.8) 



which has only one positive root given by 



Uq 



x/9ft) 2 (M 2 



;i-3ft>)4 2 



3ft)M _ y/9co 2 M 2 + (1 - 9ft) 2 )<7 2 - 3coM 
~ (l-3(0)q 2 ' 



(2.9) 
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Using this in (12.71 ) we obtain 



q 2 up-M 

c c = 75 — rr\ — ^ (2.10) 

(3(0 + l)u(f m 



for an extreme black hole solution. Since the r.h.s of ( 12.71 ) is positive, we must have «o > M/q 2 , which is 
directly derived from ( I2.9I ). The roots (w/j,wi) are such that 



u h 



< uq < u\ and M/q 2 < u\ (for c> c c ) (2.11) 



but iih may be greater or smaller than M/q 2 . 

For fixed ((0,M,q 2 /M 2 > 1) we conclude that if c < c c , the metric ( 12. 11 ) is a naked singularity and if 
c > c c , the solution is a black hole with an event horizon at r/, and an inner horizon at r\. The case c = c c 
corresponds to an extreme black hole with an event horizon at tq = \/uq with uq given by (12.9I ). 

We would formulate the last statement in terms of fixed ((0,M, c) if it were possible to solve (12.101 ) for 
q 2 /M 2 . An extreme black hole is obtained a) upon fixing ((0,M ,q 2 /M 2 > 1) and decreasing the value of 
c till c = c as done in Figured] orb) upon fixing ((0,M,c) and increasing q 2 /M 2 till q 2 /M 2 = (q 2 /M 2 ) c . 
The case b) is generally not possible because we cannot invert (12. 10b . however, for fixed c <C 1 and given 
((0,M) the extreme RN black hole in presence of quintessence is approximated by 



1 + 77^-7 + ^ 777^77 (2-12) 



(m 2 )c * 1 M 3co+1 ' " \M 3m + 

Recall that 1 >3(0 + l > 0, so that the contribution of the second term is not neglected even for relatively 
large black holes. For 3(0 + 1 = ((0 = -1/3) and < c < 1/2, it is easy to show using d24l ) or (f27Tll272l > 
that the solution is a black hole if q 2 /M 2 < 1/(1 — 2c) and that the extreme black hole corresponds to 

(w), = T^ < 0<1 - 2c<1 »' (213) 
3 Thermodynamic of RN black holes with quintessence 

In this section we only consider asymptotically flat black hole solutions which correspond to 

[(q 2 < M 2 and any c> 0) or (q 2 > M 2 and c > c c )) if - 1/3 < (O < (3.1) 



? ? 1 \/? 9 1 q 2 -M 2 



< M and - > c> or [q L > M z and - >c> 

2 / V 2 2g 2 



ifw = -l/3. (3.2) 



We first consider the case — 1/3 < CO < 0. We may take as thermodynamic state variables of these 
black holes the entropy S, the electric charge q and the parameter c. The entropy of such black holes is 
usually taken as S = nr/ 2 , which is the area of the horizon by 4; however, for simplicity of notation we 
will work with the quantity s = S/n. Using f(rn) = 0, we express the (gravitational) mass M in terms of 
(s, q, c) as follows 
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It is straightforward to show that, in presence of quintessence, M is still the (internal) energy of the black 
hole. For that purpose we need to show that (dM/ds) qc is indeed the temperature of the black hole. The 
latter is proportional to the surface gravity T = k/2 and K = d r f(r)/2. Using these equations we express 
T is terms of the entropy and mass as 



r = i 

2 



M__cf_ (3m + l)c 

7~^72 + s(3o+2)/2 



(3.4) 



Replacing M by the r.h.s of (13.31 ) reduces the r.h.s of (13.41 ) to (dM/ds) q jC . 

It is worth mentioning that the Komar mass, M^mar = lim, ^co r 2 d r f /2, diverges as 



^Komar = lim 

r— ><>=> 



M- — + 



q 1 (3fl) + l)c 



r 3c» 



However, if one subtracted the infinite background contribution due to quintessence, the remaining con- 
tribution approaches M. 

Let (Ao,Ag) be the functions 

fdM\ fdM\ fdM\ 

A ° = ( "a" ) ' A Q = ( "a" ) alon § Wlth T = ( "a" ) • 

V dq ) s,c V dc J s,q \ as J q,c 



Using (I3.3I ). we obtain 

s-q 2 + 6c(Os w - 

T = • (3 " 6) 

Using l/s = Uh 2 in (13.41 ), we see that T °< [c(3to + l)w/i 3t0 +M — ^ 2 m/,]. For the extreme black hole, c = c c 
and Uh = uq, so by (12.71 ). T = 0. 

Note that the r.h.s of (13.31) is a homogenous function of (s,q 2 ,c 2 ^ 3a)+1 ^) of degree 1 /2: 

M(As,A<7 2 ,Ac 2 ^ 3ft)+1 ' > ) = A 1/2 M(s,^r 2 ,c 2/(3(B+1) ) . 

By Euler's theorem we obtain the generalized Smarr formula 

M = 2rs + A <7 + (3a>+l)A e c (3.7) 

where Aq is the electric potential and Aq is the potential associated with quintessence. It is straightforward 
to check that changes in the thermodynamic state variables (s,q,c) by amounts (ds,dq,dc) result in 

dM = Tds + A d<7 + A e dc (3.8) 

which is the first law of black-hole thermodynamics for RN black holes in a background of quintessence. 
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For the case CO = —1/3 we use (12.41 ) where the hole has mass M' and charge q' . The horizon area 
is multiplied by 1 — 2c so that the entropy divided by n is s' = (1 — 2c){M' + a/M' 2 — q' 2 ) 2 . Let s 
v' i I - 2c) = (M' + \JM' 2 — q' 2 ) 2 , then 

M' = ^±|^ and T' = - *~' 7 '' (3.9) 
2^ 4(l-2c> 3 /2 

[where T' = (dM' / ds') q i c ]. With Aq = (dM' /dq') s > c , the generalized Smarr formula for this case takes 
the form 

M' = 2rV+Ao< ? / . (3.10) 

Particle absorption-emission: Adiabatic invariance. Without loss of generality, we assume q > 0. 
As shown in ll20l . a particle of mass m and electric charge £ > moving in the geometry described 
by (12.11 12.21) has the conserved energy for radial motion 



E=m^f(r) + (dr/dx) 2 + ^ (3.11) 

where T is the proper time. If the particle crosses the horizon r/,, this incurs changes in (M,q): dM = E, 
dq = e, which are related by (13.81 ) [here we assume that the motion of the particle does not affect the 
density of quintessence, that is we take dc = 0] 



2(dM-— ) = [l-q z u h z + 6ccou h ico+l ]dr h . (3.12) 
\ r h J 

Using /(r^) = 0, we bring it to the form 

dM - ^ = -[W-q 2 u h 2 + 3Mcou h - W+\dr h (3.13) 
rh 

where the expression inside the square parentheses is proportional to the l.h.s of (12.81) . Since, by (12.111) . 
u/i < uo, the coefficient of drp, in (13.131) is positive. Thus 

sgnfdM- n = sgnfdr/,) > (by the second law) (3.14) 
V r h ) 

which we rewrite for the charged particle 

sgn (e — — ) = sgn (dr/,) > (by the second law) . (3. 15) 

Hence, if the particle's energy is Eq = qe/r/ u the latter reaches the horizon with a zero speed: dr/dx = 
(as seen from (13.1 lb ). In this case the particle is adiabatically accreted by the hole causing no change in 
the horizon's area: dr^ = (by ( 13.151 )). 

If E > Eq the particle plunges into the hole with some kinetic energy causing the horizon to expand 
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by the amount 

(E - E )/[W+ - 3Mcou h - W_q 2 u h 2 } 

[particles with energy E < Eq cannot reach the horizon]. 

Conversely, as shown in [21'], black holes with relatively high temperature radiate electrically charged 
particles in the superradiant regime. The strong gravitational field near the horizon creates two particles 
of opposite charges. The total energy lost by the field is greater than the energy carried by each particles 
and in any casqj dM < 0. If dq > 0, that is the particle with positive charge falls into the hole and the 
particle with negative charge escapes to infinity, the l.h.s of (13.141 ) is negative (since dM < 0) leading to 
dr/, < which is not favored by the second law. Thus, by the second law of thermodynamics, the black 
hole may radiate charges of only the same sign as its own charge q reducing its charge upon receiving 
electric charges of opposite sign to its charge so that dq < 0. The process may continue till the entropy 
reaches its maximum values or proceed adiabatically (i.e. the charge e of the falling particle is such that 
dM = qe/rf, < 0) till the temperature drops. 

Adiabatic charging of an RN black hole. We consider an ordinary RN black hole (BHn), where c = 0, 
along with another RN black hole surrounded by quintessence (BHq) for given and fixed values of {00, c) 
with — 1/3 < CO < 0. Each hole has mass Mo and charge qo satisfying Mo/qo > 1 as in Figured] (a). The 
horizons, r + of (BHn) and of (BHq), are such that «/, < u+. 

Let r e denotes r + or r/, (u e = \/r e denotes h + = l/r + or «/, = lA/ 7 ). If the two black holes are 
charged adiabatically, r e remains constant (throughout the rest of this section u + and w/, are then taken as 
constants), so that by (13.121 ). dM = qdq/r e leading to 

M=^-q 2 + D (D = M () -y«?o 2 ) (3-16) 

where the constant D > since u e < M /q 2 as shown in Figure \T\ (a). Using ( 12.51 ) it is easy to show 

For (BHn): 2Du+ = 1 (2D = r+) (3.17) 
For (BHq): 2Du h = 1 - 2cu h 3m+l < 1 . (3.18) 

While the two holes are being charged (at the same rate), the point M/q 2 (Figure [J (a)) moves pro- 
gressively to the left since 

M u e D 
q 2 2 q 2 

is obviously a decreasing function of q 2 . As M/q 2 meets first u + , (BHn) turns into an extreme black hole 
and the charging process ends [since (13.141 ) is no longer valid] by cumulating the total charge and mass 

For (BHn): u + = ^± + -2r =¥> q, = r + , M, = r+ (3.19) 
2 q t l 

2 Excluded is the case where both particles fall into the hole. If we neglect radiations par the particles, this corresponds to 
dM = 0, dq = and dr h = 0. 
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as expected (the subscript "t" for total). What happens at this moment to (BHq)? Let us look at the 
derivative of M/q = u e q/2 + D/q: 

d(M/q) _u e D 
dq 2 q 2 

This is zero if q 2 = 2D/u e . For (BHn), using (13.171 ) along with u e = u + , we obtain q = r + = q t , which 
means that M/q reaches its minimum value, 1, at the end of the charging process. For (BHq) the situation 
is quit different: At the moment M/q 2 meets u + , the function M/q is still decreasing, so that (BHq) 
cumulates more charge and mass than (BHn), and reaches its minimum value at the moment M/q 2 meets 
Uh with 

For (BHq): (M/q) min = y/2Du h < 1. (3.20) 

As the point M/q 2 passes uh, the ratio M/q starts to increase but remains smaller than 1. The charging 
process ends when (BHq) turns into an extreme black hole at the moment the point uo (Figure [U(c)) meets 
The total charge and mass of (BHq) are then given by (use (12.81 ) and (13.161 )) 

(1 - 3co)q t 2 u h 2 + 6M,cou h - (3(0 + 1) = and 2M, = u h q t 2 + 2D 

leading to 

For (BHq): q t = r h ^\ + 3(o(\-Du h ) < r h (3.21) 
For (BHq): 2M t = r h [(l - 3(o)2Du h + (3(0 + 1)] < 2r h . (3.22) 

It is easy to check that: 

For (BHq): M t /q t < 1 , r + < q t < r h and r + <M t <r h . (3.23) 
4 Thermodynamic local stability 

In this section we only consider thermodynamic processes for which c is held constant. The second law 
of (black holes) thermodynamics governs all the criteria for stability. These criteria, in form, depend 
on how we consider the system (here the black hole). We shall consider canonical and microcanonical 
ensembles (CE, ME, respectively). The CE will consist of the black hole in equilibrium with its thermal 
radiation, treated as a reservoir (heat bath) at constant temperature and the ME will be the case of a black 
hole isolated from its surroundings. It has been shown in many applications that the thermodynamic local 
stability of CEs may equally be treated by classical thermodynamic approaches based on the Hessian 
matrix of the entropy or, equivalently, of the energy [4 -6] or by the TPM [8 ] derived by Poincare Q. 
Due to non-additivity of entropy and mass in general relativity (9], the classical thermodynamic approach 
does not apply to isolated black holes, for which we shall then apply the TPM. 
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4.1 CE: Classical thermodynamic approach 

We assume that the hole is immersed in a thermal bath at constant temperature T. Applying the classical 
thermodynamic approach, we denote by ds b the change in the entropy of the bath. Then any possible 
change in the state of the system requires ds + dsb > 0. Conversely, if fluctuations try to take the system 
out of equilibrium with the reservoir, that is if 

ds + ds b <0 (4.1) 

for allowed changes in the system's and reservoir's parameters, the system cannot leave the current state 
which is said to be in stable equilibrium with the bath. The inequality (14.11) is the condition from which all 
criteria for local stability of systems in contact with reservoirs are derived. Since any process is considered 
reversible for the (huge) reservoir, from (13.81 ) it follows that dsb = (dMb —Aodqb)/T (dc = 0). Exchanges 
between the system and the reservoir obey the conservation rules: dMb = —dM and dq b = —dq. Using 
these equations in (14.11) leads to 

dM >Tds+A dq. (4.2) 



In ( 14.21) we are using the mass-energy as a fundamental thermodynamic quantity, instead of the entropy, 
because it is not possible to reverse (13.31 ) and express s in terms of M. 

By the first law, dM = Tds + A(,dq, the condition ( 14.21 ) is not sensitive to first order changes in the 
allowed parameters. Keeping up to second partial derivatives of M with respect to the extensive param- 
eters (s,q) in the Taylor expansion of M, we obtain the equivalent condition (which may be written as a 
Hessian matrix) 

l(d 2 M\ l2 (d 2 M\ , , \(d 2 M\ l2 n ,„ „ N 

where all first order terms canceled out. With 



d 2 M\ _T /d 2 M\ _(dT\ /d 2 M\ _ /dA 
ds 2 ) q,c C ? ' \dsdqJc \dq)s,c' 1 \ dq 2 ) s,c V dq 



(4.4) 



the condition (14.31) implie; 



(T and Aq have been defined in the equation preceding ( 13.51) and their expressions are given in ( 13.51 13.71 )), 



(4.5) 



/dA \ o an d( dT Y T ( dA ° 

' V dq J s,c \dq)s,c C q \ dq 



3 We would like to show the analogy with classical thermodynamics via the following correspondences where (P , V) are the 
pressure and volume of the classical system: dg — ¥ —dV, Aq — > P and C q —¥ Cv ■ In a similar way, (dq/dAo) s , c /q, which may 
be called the factor of adiabatic charge, corresponds to the adiabatic compressibility —(dV/dP)s/V. Similar terminology has 
been used in |22|. 
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Here C q is the specific heat at constant charge 



/dM\ ,/d 2 M 



q ~\dTJ q ,c~ '\ds 2 sq,c 
2s(s-q 2 + 6c(Os w -) 
' 3q 2 -s-6cco{2 + 3(o)s w - ' 



(4.6) 



Note that if T is constant (the case where the hole is immersed in a reservoir), the l.h.s of the third 
inequality in ( 14.51 ) is zero. 

The mass M, given by (13.31 ). is supposed to be positive, so we have the extra condition 

2j~sM = s + q 2 - 2cs w - > . (4.7) 

Usually, in classical thermodynamics, all relevant thermodynamic quantities pertaining to the reser- 
voir (Mb,qt,,- ■ ■), but the temperature, are allowed to fluctuate. In this paper we shall allow T to fluctuate 
too and investigate separately the cases T constant and T fluctuating. 

4.1.1 T constant 



Condition (14.31 ) or (14.51 ) reduces to C q > and T > 0. Since s > and C q °< T, we have to solve simulta- 
neously 

(a): s - q 2 + 6ccos w - > and (b): 3q 2 - s - 6cco(2 + 3co)s w - > . (4.8) 
For c = we recover the conditions for ordinary RN black holes: 

q 2 < s < 3q 2 (4.9) 

with s = (M+ \JM 2 — q 2 ) 2 leads to the known conditions for local stability |6] 

3/4 <q 2 /M 2 < I. (4.10) 

Now back to quintessence case c > 0. For — 1/3 < ft) < we have 1/2 < W_ < 1 and 1 < 2 + 3w < 2. 
In order to solve the inequality (14.81) (a), we consider the line y(s) = s — q 2 and the concave-down curve 
y(s) = —6c(Os w - (the graph of which is similar to that ofy = y/s), it is easy to see that they do intersect 
at one and only one point si (ft)) > q 2 . Thus T > if s > s\ (ft)), which solves (14.81 ) (a). Similarly the line 
y(s) = s — 3q 2 and the concave-down curve y(s) = —6cco(2 + 3(o)s w - intersect at one and only one point 
S2(co) > 3q 2 . Thus (14.81 ) (b) is solved by s < ^(w). The black hole is locally stable if 

si (to) < s < s 2 ((o) and M > 0. (4.11) 

The first condition in ( 14.111 ) is a generalization of ( 14.91 ) where q 2 and 3q 2 have been shifted to the right 
(Ji(co) > q 2 , S2{(o) > 3q 2 ). The problem, however, cannot be solved as in (14.101) since {s\,s,S2) are not 
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known explicitly in terms of (M,q, co,c). 

In the case of ordinary RN black holes (c = 0), the mass, 2^/sM = (s + q 1 ), is a positive function of 
the entropy. In the case of RN black holes surrounded by quintessence, this is not guaranteed a priori and 
we need the solve the extra inequality ( I4.7I ). Consider the line y{s) = s + q 2 and the concave-down curve 
y{s) = 2cs w - . For fixed (ft), c), if we choose q 2 > q 2 = W+ (2cW- [for ordinary RN black holes 

(c = 0), q 2 = 0], the line and the curve do not intersect and M > for all s and if q 2 < q 2 , they meet at 
two points 54(a)) < 53(a)) where M > for s > 53(a)) or s < 54(a)). 

Figure[2](a) shows plots of the curves 5i = 5i (ft)) (red), s 2 = s 2 (oo) (blue) and 53 = 53(a)) (green) for 
— 1/3 < ft) < 0. We have chosen q = c = 1 in such a way that 54(a)) < 5i (ft)) so that there is no need to 
plot the curve 54 = 54(a)). T > for s > S\, C q > for s\ < s < s 2 and M > for 5 > S3. Three isotherm 
curves are shown in Figure 12(a): T = (red), T = l/(8ff\/2) (brown) and T = 9/(47il0 3 / 2 ) (cyan). The 
plane region of local stability, identified by T > 0, C q > and M > 0, is the region enclosed by the green 
curve, blue one and the s-axis for a>o < CO < where (a>o,so) = (—0.19,21.37) is the intersection point of 
the curves s 2 = 52(a)) (blue) and S3 = S3 (a)) (green). For g held constant, (a^^o) lies on the curve given 
by 

1- a)(3a) + 2) 2 

5 = 3- — rr<7 (Vc). 

l + 3a)(3a) + 2)^ v ; 

The region where the hole is unstable is above the blue curve for 00q < 00 < 0. Since C q diverges on the 
blue curve, this latter determines a limit for second order phase transition. 

The conditions in (fl~8l) and (l4~7l) are rewritten as T > : q 2 < s + 6c(0s w -, C q > : g 2 > (j/3) + 
2ca)(2 + 3oo)s w - (with T > 0) and M > : > — s + 2cj , respectively. Figure [2] (b) shows plots of 
the surfaces q 2 = s + 6coos w - (red), q 2 = (s/3) +2c(t)(2 + 3(o)s w - (blue) and q 2 = —s + 2cs w - (green) 
for c = 1. The curves s\ = s\(oo) (red), s 2 = ^(a)) (blue) and S3 = 53(a)) (green) are projections of 
intersections of these surfaces with the plane q 2 = 1 . The physical space region is all of the space region 
bounded above by the red surface and below by the green one if there q 2 > or the subregion of it where 
q 2 > 0. If the point p representing the thermodynamic state of the black hole lies between the red and the 
blue surfaces and above the green one (M > 0), the hole is thermodynamically stable against fluctuations 
in (s,q): fluctuations are "entropically" suppressed. Since T does not fluctuate, the stability condition is 
subject to the further constraint: r.h.s of ( 13.61 ) = constant. For fixed c, this defines a new surface in the 
space (oo,s,q 2 ) of parameters. Stability concerns only those states p of the hole which lie on the segment 
of this new surface T = constant which is sandwiched by the red and blue surfaces. The space region of 
instability is bounded above by the part of the blue surface where q 2 > and below either by the green 
surface if there q 2 > or by the plane q 2 = 0. 

For 00 = — l/3(0<c<l/2)we have seen that the solution is an ordinary RN black hole, given 
by (12.4b . with new mass M' and charge q'. Using ( 13.9b we obtain 



[Recall that in this case s = s' /( 1 — 2c) = (M' + \JM' 2 — q' 2 ) 2 ]. Thus the conditions for local stability are 
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(Si, S 2 , S 3 ) [til) 

30 r 




(a) 




a (b) 



Figure 2: (a) Plots of Si = s\(co) (red), s 2 = s 2 {co) (blue) ands 3 = s 3 (co) (green) for - 1/3 < CO < and q = c = 1. r > for 
i - > ji, C ? > for s\ < s < S2 and M > for j' > .53. The red, brown and cyan curves are the isotherms T = 0, T = 1/ (8;r\/2) 
and T = 9/(4s:10 3 / 2 ). The curves 52 = S2(a>) (blue) and J3 = .53(0)) (green) intersect at the point (oo,^o) = (—0.19,21.37). 
The black hole surrounded by quintessence is locally stable in the region between the green and blue curves for C0q < CO < and 
unstable in the region above the blue curve for C0q < CO < 0. Limited to these values of CO, the blue curve determines the points 
where C q diverges, there the black hole undergoes a second order phase transition, (b) Plots of the surfaces q 2 = s + 6ccos w - 
(red), q 2 = (s/3) + 2cft)(2 + 3co)s w - (blue) and q 2 = — s + 2cs w - (green) for c = 1. T > below the red surface, C q > below 
the red surface and above the blue one and M > above the green one. 

the same as in (14.91 ) but with M' and q' , leading to 

3/4 < (\-2c)q 2 /M 2 < 1. (4.13) 



By (13.21 ) we see that black holes with either q 2 < M or q 2 > M 2 are locally stable provided condi- 
tions ( 14.131 ) are satisfied. Black holes with 3M 2 /4 > (1 - 2c)q 2 are unstable for C' q < 0. In the (M 2 ,q 2 )- 
plane, the line M 2 = 4(1 — 2c)q 2 /3 determines a phase transition of second order where C' diverges. 



4.1.2 T fluctuating 

In processes where T fluctuates, the conditions in (14.51) need be solved simultaneously. The second 
condition in (14.51 ) is already satisfied; using (13.51 13.61 14.61) . the third one leads to 



q 2 -s-6c(0{2 + 3(0)s w - 
3q 2 -s-6ca>(2 + 3(0)s w - 



>0. 



(4.14) 



Note that if the numerator of the fraction in ( 14.141 ) is positive, the denominator, which also appears in the 
expression of C q , too will be positive. Thus, all we have to solve are the conditions (14.7114.81 (a)) and the 
new one q 2 — s — 6c(o{2 + 3d))s w - > 0, all grouped into one system, respectively 



s + q 2 >2cs w - (M>0) 
s-q 2 > -6cC0s w - (T >0) 
s-q 2 < -6cw(2 + 3co)s w - . 



(4.15) 
(4.16) 
(4.17) 
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We first consider an ordinary RN black hole (c = 0). The last two conditions (14.161 14.171) lead to 
s = q 2 , then, by the first one, s = q 2 > 0. With s = (M + \JM 2 — q 2 ) 2 , we obtain M = q. Thus, only 
the extreme RN black hole is stable against fluctuations in (s,q,T). The first law takes the simple form 
dM = dq: for any change in q there corresponds an equal change in M so that the hole remains extreme 
for any dq. 

In the generic case c > 0, -l/3<O)<0we have 1 < 2 + 3ft) < 2. The solution of (14.171 ) is similar 
to that of (14.161) : the line y(s) = s — q 2 intersects the concave-down curve y(s) = — 6cft)(2 + 3co)s w - at 
one and only one point S5(a>) > q 2 [note that ss(co) < ^(ft*)]- Now, because 2 + 3(0 > 1, the graph of 
y(s) = — 6cft)(2 + 3(o)s w - is above that of y(s) = —6ccos w -, the intersection of which with y(s) = s — q 2 
determines the point s\(co) > q 2 as discussed earlier. This leads to si(ft)) < s 5 (ft)). Thus, the three 
conditions (14. 15l to 14.171) are solved by 

Ji(ft)) < s< s 5 (cd) and M>0. (4.18) 

Allowing fluctuations in T, the condition for local stability has been narrowed since ss((o) < S2((o) 
(compare with (14.1 II) ). This is the region bounded below by the red curve and above by the magenta 
one of Figure [3] (a), which shows plots of the curves s\ = si(ft)) (red), S2 = *2( ffl ) (blue) and s 5 = s 5 (co) 
(magenta) for— 1/3 < ft) < 0, c = 1 and q = 5. The three curves do not intersect and admit the line 
ft) = —1/3 as vertical asymptote. For fixed ((0,q), as the entropy increases from the red line which 
defines the states of extreme black holes, the state of the hole crosses the magenta line and becomes 
unstable regarding fluctuations in T. However, this transition cannot be qualified a first order phase 
transition since the entropy is continuous there (no jump in s) and the phase of the hole is the same; rather 
one might qualify it a "behavioral" change or transition. In the region bounded below by the magenta 
curve and above by the blue one, the hole is, however, stable regarding changes in s, q, as we have seen 
in the previous subsection. So this is a phase transition from thermodynamic states, which are stable 
against fluctuations in (s,q,T), to states which are stable against fluctuations in (s,q) only. If the entropy 
continues to increase, the hole undergoes the above-mentioned second order phase transition by crossing 
the blue line where C q diverges and changes the sign. 

Figures[3](b,c) show plots of the surfaces q 2 = s + 6ccos w - (red), q 2 = s + 6c(o(2 + 3co)s w - (magenta), 
q 2 = — s + 2cs w - (green) and the plane q 2 = 25 (yellow) for c = 1 [for clarity, the blue surface q 2 = 
(s/3) + 2cft)(2 + 3co)s w - is not shown]. The plane q 2 = 25 intersects the red surface, which represents 
the states of extreme black holes, along a curve the projection of which on the (ft), s)-plane is the red curve 
of Figure [3] (a). The other curves in Figure [3] (a) are also projections of intersections of the plane with 
the corresponding surfaces. For fixed (ft), s), as q 2 decreases, along a vertical line (Figure [3] (c)), from 
its value on the red surface to its value on the magenta one, the hole remains stable against fluctuations 
in T till the line crosses the magenta surface and the hole undergoes a behavioral change. As the charge 
decreases again, the state of the hole crosses the blue surface (not shown) and undergoes a second phase 
transition of second order. 

We have thus shown that black holes with low entropy, or high charge or both are stable against fluc- 
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Figure 3: (a) Plots of s\ = s\(co) (red), S2 = S2(®) (blue) and s 5 = s 5 (co) (magenta) for -1/3 < (0 < 0, c = 1 and q = 5. (b) 
The space region of stability against fluctuations in both (s, q, T) is above the magenta surface and (c) below the red one and 
where the green surface is below both of them. 

tuations in (s,q, T), and black holes with relatively high entropy, or low charge or both are stable against 
fluctuations in (s,q) only. There is a behavioral change between the states of these holes maintaining the 
sign of C q > 0. Finally, black holes with high entropy, or low charge or both are unstable against allowed 
fluctuations and have C q < 0. 

For (o = -1/3 (0 < c < 1 /2), using Eqs. ( 13.11114.121) we bring the third condition in (|43T ) to 

a' 2 -s 

->0. (4.19) 



3q' 2 

The problem is similar to that of ordinary RN black holes with M and q replaced by M' and q'\ Solv- 
ing (14.191 ) together with T' > and C' q > we obtain M' = q' or M = a/1 — Icq. Only the extreme black 
hole is stable against fluctuations in (s,q). The first law takes the simple form AM = x/T— 2cdq: for any 
change in q there corresponds a proportional change in M so that the hole remains extreme for any dq. 
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4.2 ME and CE by the TPM 

We apply the TPM to the case of an isolated RN black hole surrounded by quintessence, which is the ME 
and extend the analysis to the case where the hole is immersed in a heat bath, which is the CE treated 
in the preceding subsection by the classical thermodynamic approach. The aim of applying the TPM to 
the CE is to rederive in an elegant way the results of subsubsections 14.1.11 and 14.1.21 using the powerful 
method of Poincare Q. 

In the TPM one uses the following Massieu functions *P: (s, —F/T, —G/T, • • •) or, equivalently, the 
thermodynamic potentials x ¥: (M,F,G, • • •) |[8] depending on the ensemble, where F = M —Ts and G = 
F — A$q are the Helmholtz and Gibbs free energies, respectively. For our ensembles the corresponding 
potentials at equilibrium are: 

ME: V =M(s,q) (4.20) 
CE: *P = F(T,q) if T constant (4.21) 
CE: ¥ = G(T,A Q ) if T fluctuates . (4.22) 

The TPM consists in plotting the planar curves d^/dx (all other variables kept constant) against x, 
which is some control parameter. The variable d^/dx is called the conjugate of x with respect to x ¥. 
These curves are called linear series of equilibrium. Changes of equilibrium (from stable to less stable 
to unstable and conversely) occur at points where the curves have vertical tangents or bifurcations. If the 
curve has no vertical tangents, as in Figure @] (a), then all points on the curve have the same degree of 
stability: If it is known that a point on the curve represents a stable equilibrium state, then all the points 
on the curve represent similar states. If the curve, with vertical tangent, is concave left near the vertical 
tangent (as shown in Figure [4] (b,c)), then all the points on the branch of the curve where the slope is 
negative near the vertical tangent (upper branch in Figure @] (b,cV) are more stable than the points on the 
other branch where the slope is positive near the vertical tangentj (lower branch in Figure H](b,c)). 

For the ME, the control parameter is the entropy s. Since (d x ¥/ds) q = (dM/ds) q = T we have plotted 
in Figure|4](a) the series of equilibria T(s), with q constant, for both the isolated RN black hole surrounded 
by quintessence (fitted line) and Schwarzschild black hole (dotted line). The curves approach each other 
as s — > oo. Knowing that the isolated Schwarzschild black hole is locally stable JH, we conclude that the 
isolated RN black hole surrounded by quintessence is at least stable for large values of s. But since there 
is no change of stability on the series of equilibrium (no vertical tangents on the curve y = T(s)), the hole 
is then stable for all s. 

The case of the CE has been split into tow subcase: T constant (subsubsection 14. 1 7Tb and T fluctuating 
(subsubsection |4TL2l ). In the former case, *P = F(T,q) and T is a control parameter. Using the first law, 

4 The slope need not be of the same sign along a given branch. The rule is valid if is one of the potentials (M,F,G- ■ ■ ), 
which are minimum for stable equilibria. For *F chosen from the list of Massieu functions (s, —F/T.—G/T, ■■■ ), which are 
maximum for stable equilibria, then if the curve is concave left, all the points on on the branch of the curve where the slope is 
positive near the vertical tangent are more stable than the points on the other branch where the slope is negative near the vertical 
tangent. 
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Figure 4: (a) (q = l,c = I, CO = -0.1): Plot of T vs. s at constant q for both the isolated (ME) RN black hole surrounded by 
quintessence (fitted line) and Schwarzschild black hole (dotted line), (b) (q = 1 , c = 1 , CO = — 0. 1 ) : Plot of —5 vs. T at constant 
q for the RN black hole surrounded by quintessence immersed in a heat bath (CE). (c) (T = 0.05, c = 1,0)= —0.1): Plot of Ao 
vs. q at constant T for the RN black hole surrounded by quintessence immersed in a heat bath (CE). It is a parametric plot for 
si < s < s 2 , s 2 <s< s , sj = 0.329134, s 2 = 7.191, s = 14.65, q = V 's-ATs^l 1 + 6ccos w - and A = q/^/s. 

dM = Tds +Arjdq, we obtain 



Hence, (d x ¥/dT) q = (dF/dT) q = —s. We need to plot the curve — s(T) for q constant. This is the same 
curve plotted in Figure |4] (a) but rotated 90° clockwise as shown in Figure @](b). Notice that the point A 
where —s(T) has a vertical tangent is the same point where T(s) has a horizontal tangent. Using (13.61) 



we see that (dT /ds)q jC = at the points S2 where C q diverges. Thus we reach the same conclusion as in 
subsubsection 14. 1 . Tl For fixed (q,c,co), the stability breaks at the point A where there is a vertical tangent 
(where C q = °°), the points where T < T2 = T(s2) and s\ < s < S2 are stable (negative slopes) and the 
points where T < T2 and s > S2 are unstable (positive slopes). 

For the ordinary RN black hole (c = 0), si = q 2 , S2 = 3q 2 and T% = \j (6\/3^)- This is the maximum 
temperature beyond which the hole cannot resist evaporation. 

Using (14.231 ). we have (d x ¥/dq)x = (dF /dq)j = Aq. So we could also plot Ary(q) for T constant, 
using thus q as a control parameter. This is a parametric plot, shown in Figure 0](c), parameterized by s. 
The entropy increases along the upper branch (stable) from si at O to S2 at B, it continues to increase along 
the lower branch (unstable) from B back to O. To show that, at B, s = S2, which is the point where C q 
diverges, we rewrite (13.61) as q = Vs-4Ts 3 / 2 + 6ccos w -. With T constant ((c, ft)) are also kept constant), 
at B we have dq/ds = leading to s — 6Ts 3 ^ 2 + 6c(dWs w - = 0. Eliminating T by (13.61) . the remaining 
equation is again 3q 2 — s — 6cft)(2 + 3(o)s w - = (compare with (I4.24I )). the solution of which is s = S2- 
Here again we reach the same conclusion as in subsubsection 14.1.11 

If, in the CE, we assume that T fluctuates, the appropriate potential is the Gibbs function *P = G(T, Ao) 
leading dG = —sdT — qdA$. Using Ao as a control, we obtain {dm / dA^jj = (dG/dArj)T = —q. So we 
need to plot the curve —q(Ao) for T constant. This is the same curve plotted in Figure [4] (c) but rotated 
90° clockwise. So the horizontal tangent at C becomes a vertical one and the (upper) branch from si 
at O to S5 at C is the stable one if fluctuations in T are taken into consideration [negative slopes near C 



dF = -sdT+A Q dq. 



(4.23) 




(4.24) 
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after revolving the curve in Figure @](c) 90° clockwise]. The branch from C through B back to O is the 
unstable one regarding fluctuations in T. So we just re-derived the conclusion made in subsubsection |4. 1.2] 
(compare with (14.1 81 )). All we need to show is s = s$ at C. At C, dAo/ds = which is the same as 
2Ts i l 2 + 3cft)( 1 + 3co)s w - = 0. Eliminating T by ( 13.61 ), the remaining equation is again q 2 — s — 6cft)(2 + 
3co)s w - = (compare with (14.171) ). the solution of which is s = s$. 

5 Phase transition via geometric methods 

In this section we briefly present the results of thermodynamic stability as derived by the two geometric 
approaches: 1) Geometrothermodynamics (GTD) lfl5l and 2) Liu-Lu-Luo-Shao (LLLS) method ifTBl . 

5.1 Geometrothermodynamics 

Let *P be a thermodynamic potential and (E a ,I a ) the set of associated extensive and intensive variables. 
We define a (2n + 1) -dimensional space T whose coordinates is the set Z A = {m,E a ,I a } where A : — > 
In and a : 1 — > n. Together with the Gibbs 1-form = d*F — 8 a bI a dE h , (T,0) make up the (2n + 1)- 
dimensional contact manifold of metric G AB (Z C ), which is the thermodynamic phase space [15]. Here 
we assume that T is differentiable and that satisfies the condition A (d©)" ^ 0. The subspace E C T 
of equilibrium states is defined by the map: <p : E — > T such that the pullback vanishes at 0: <p*(0) = 0. 

The singularities of the curvature scalar Rgtd of E determine the points or states where there are 
second order phase transitions of the thermodynamic system. The metric of the space E is given by ll23l 



where r] ac i = (— 1,1,..., 1). This metric is invariant under Legendre transformations lfl5l . 

We consider an RN black hole surrounded by quintessence and immersed in a heat bath at fixed 
temperature, we assume no fluctuations in c (dc = 0): This is the CE subject to (14.3I ). It is then convenient 
to use the mass ( 13.31 ) as the thermodynamic potential. We can write the matrix d5.ll ) in a metric form 

„, / dM dM\ ( d 2 M , d 2 M - 

{[s-3q 2 + 6cco{2 + 3(o)s w -}ds 2 + 8s 2 dq 2 } (5.2) 



8s 5 /2 

and the associated curvature scalar reads 



X 



^GTD--^ 



{-3q 2 + s - 18cW 2 ^-) 2 + \Sq 2 [-3q 2 +s + 6c(o{2 + 3w)s w -^ 



+ 



64s 3 

X 2 [2s + 3c(o(\0 + 2\(O + 9(O 2 )s w -] 3 r q 2 - s + 6c(0 2 (5 + 6a)s w - 



2Y I 16s 9 / 2 



+ 2F 2 r 16s 5 /2 



X 2 (V ~ j + 1 8c(Q 2 s w - )[-5g 2 + s + 2cco(S + 1SCQ + 9(Q 2 )s w - } 
+ 128? 



} (5-3) 
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where 



C a 8s 5 / 2 

X = — = —a ; r^r (5-4) 

T 3q 2 -s-6cco(2 + 3(o)s w - 

3q 2 + s + 6ccos w - 

Y=sT+qA = q ^ /a- • < 5 " 5 ) 

It is obvious that Rqtd diverges at the points where X diverges corresponding to C q = °°. A divergence 
in the value of C q announces a second order phase transition as derived in subsection l4.1.1l 7?gtd diverges 
at Y = too but this equation, based on the analysis made in subsections 14.1.11 and 14.1.21 represents no 
physical effect. This is a pathologic effect due to the fact that the metric (15.21 ) becomes singular at Y = 0. 



5.2 Liu-Lu-Luo-Shao method 

The second geometric method was developed later by Liu-Lu-Luo-Shao (LLLS) [18]. The method relies 
on the same principal as the first one in that any singularity in the curvature scalar of the associated metric 
signals a second order phase transition of the system. The metric of this space is simply the Hessian matrix 
of the Helmholtz free energy which, in the CE, is given by 

,2 _ »r A . lAl _ A „_ BT AAi (dM E^^.^j 



d/j t r s = -&T&s + dA d<3 = — —ds z + - — dsdq + -^-dq z . (5.6) 

ds \ as oq J aq 

Using the expressions ( 13.51 ) and (13.61 ), ( 15.61 ) becomes 

d/ LLLS = ^{[* " 3tf 2 + 6cco(2 + 3co)s w -]ds 2 + 8, 2 d^ 2 } (d/ 2 TD = yd/ 2 LLS ) . (5.7) 

This metric provides the following curvature scalar 

[2s + 3c(o(l0 + 2\(O + 9(O 2 )s w -} v2 

«LLLS = 22^4 ' ' ' 



Here again we clearly see that the divergence of the curvature scalar ( 15.8b corresponds to that of C q , thus 
by the LLLS method we are too able to locate the points where the second order phase transition takes 
place. 



6 Conclusion 

We have seen that asymptotically flat RN black holes surrounded by quintessence have higher entropies 
than ordinary ones. For a given value of the quintessence density c, they commulate more electric charges 
than ordinary RN holes before they become naked singularities. For c <C 1 , the maximum relative cumu- 
lated charge (q 2 —M 2 ) /M 2 is proportional to c/M 3ffl+1 , if -1/3 < 03 < 0, or to 2c/(l - 2c), if 03 = - 1 /3. 

Taking c as a thermodynamic variable, as some works do for the cosmological constant H241I251 . we 
have obtained the generalized Smarr formula and the first law of thermodynamics. 
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As one charges adiabatically an ordinary RN black hole, it cumulates mass and charge till their total 
values equate the radius of the horizon, which remains constant during the process. If the RN black hole 
is surrounded by quintessence then it will cumulate more mass and charge, however, with their totals 
never exceeding the radius of the horizon. 

Applying the classical thermodynamic method and restricting ourselves to the CE we have obtained 
generalized conditions for local stability of RN black holes surrounded by quintessence. These conditions 
are the shifting of the known ones for ordinary RN black holes if fluctuations in T are not allowed; If T 
fluctuates, the same conditions apply with their upper limits constrained to lower values. We have reached 
the conclusion that , allowing fluctuations in (s,q,T), only black holes with low entropy, or high charge 
or both are stable, while black holes with relatively high entropy, or low charge or both are stable against 
fluctuations in (s,q) only. Between these and the other black hole states there is a behavioral change 
maintaining the sign of C q . 

We have completed the analysis of stability by applying the TPM and obtained an upper limit for the 
temperature (T2 = Tfa)) beyond which the CE is no longer stable, if fluctuations in T are not allowed. 
For ordinary RN black holes T2 = l/(6\/3#). If T fluctuates, we have obtained an upper limit for the 
electric potential on the horizon (Ao(js)) beyond which the CE is no longer stable. Another general result 
we could derive is that isolated black holes, which constitute MEs, are stable. 

By the two geometrical method (GTD, LLLS) we could also determine the states corresponding to a 
second order phase tr ansition. 

The stability of non-asymptotically flat solutions is more involved and constitutes the matter of a 
subsequent work along with the case where c fluctuates. 
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